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Abstract 



CH ■ It is shown that on the de Sitter spacetime global quantum modes 

can be prepared by a global apparatus which can measure QED tran- 

,—4- | sitions as in special relativity. The difference is that the gravity 

\£> ■ eliminates the mass-shell constraints giving rise to QED transitions 

with non-vanishing amplitudes even in the first order of perturba- 
tions. These amplitudes vanish in the flat limit where the mass-shell 

f— s ■ constraints are restored. Of a special interest could be the first order 

amplitudes of a new effect of dynamical particle creation due to the 
expansion of the de Sitter universe. 
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1 Introduction 

The theory of quantum fields on curved spacetimes is faced with the sensitive 
problem of separating the subspaces of the particle and antiparticle quantum 
modes. In special relativity the Poincare isometries leave these subspaces in- 
variant determining thus a stable vacuum state. In general relativity the 
curved spacetimes could have no isometries such that the criterion of sepa- 
rating the particle and antiparticle modes may be local. Then the vacuum 
becomes unstable in the sense that the quantum modes of two different frames 
may be related through Bogoliubov transformations which mix the particles 
and antiparticles modes among themselves |Tj. This is the mechanism of cos- 
mological particle creation giving rise to local effects in accelerated frames as, 
for example, the Unruh [2J [3] and Gibbons-Hawking [I] ones which may be 
observed using the local particle detectors proposed by Urunh [2] and DeWitt 
[5]. However, in the special case when the curved manifold has isometries 
able to take over the role of the Poincare symmetry we can define a global 
apparatus which prepares and measures simultaneously the quantum states 
as in the flat case, without supplemental local detectors [6]. Then the vac- 
uum state prepared by this apparatus is stable and the isometry generators 
become the principal operators (or observables) of the quantum theory. 

The theory of the fields with spin on local- Minkowskian curved manifolds 
has a specific structure since the spin half can be correctly defined only in or- 
thogonal local (non-holonomic) frames [3 [8]. Therefore, the Lagrangian the- 
ory of the matter fields must be written in local frames assuming that this is 
tetrad-gauge covariant. In this framework the gauge group (of the Minkowski 
metric) L\ C SO(l,3) and its universal covering group, SX(2, C), become 
crucial since their (non-unitary) finite-dimensional representations induce the 
covariant representations according to which the matter fields transform un- 
der isometries [9]. The generators of these representations are the differential 
operators produced by the Killing vectors associated to isometries according 
to the generalized Carter and McLenagan formula [TQl 13 E] ■ Thus in curved 
spacetimes with isometries there are operators globally defined on the entire 
manifold, independent on the natural or local frames one uses. Among these 
operators that commute with those of the field equations one can select the 
sets of commuting operators able to define the quantum modes as common 
systems of eigenfunctions. With this method the quantum modes are glob- 
ally defined on the entire manifold and the vacuum is stable being invariant 
under isometries. Then the quantization can be done in a canonical manner 



like in special relativity as we have shown in Refs. [12]- [I5]- We stress that 
our conjecture is closer to the genuine principles of quantum mechanics but 
this does not exclude or contradict the cosmological particle creation which 
is observed in accelerated frames using local detectors [31 [16] . 

Our method is helpful on the de Sitter spacetime where all the free field 
equations can be analytically solved while the 5*0(1,4) isometries offer us 
a large collection of operators commuting with the operators of these equa- 
tions. A special problem of this geometry is that the only time-like Killing 
vector is not time-like everywhere. This generates some doubts concern- 
ing the possibility of defining correctly the energy operator, the in and out 
fields and the scattering operator, S [17J. However, we have shown that 
there are no real impediments since the time-like Killing vector keeps this 
property everywhere inside the light-cone where an observer can perform 
physical measurements |18j . We remind the reader that in this geometry the 
light-cone domain of any observer remains behind his horizon such that the 
measurements can be done as in the flat case without to be affected by the 
presence of the observer's horizon. Thus quantum transitions can be mea- 
sured and their transition amplitudes can be calculated using perturbations 
in terms of free quantum modes. Since all the basis-generators of the so(l, 4) 
algebra have a well-defined physical meaning [18] , we can select various sets 
of commuting operators determining different sets of free quantum modes. A 
particular feature of this algebra is that the energy and momentum operators 
do not commute to each other and, consequently, the dispersion relation is 
absent and there are no mass-shells [6j [18] . This affects the measurements of 
the energy and momentum which are diagonal in different bases, called here 
the momentum and respectively energy representations. Under such circum- 
stances, we derived the principal quantum modes of the free Dirac [121 IT3] . 
Proca [E] and Maxwell [T5] fields for which we used the same method of 
canonical quantization with a stable vacuum state of the Bounch-Davies type 
[j~9] . Moreover we obtained the expressions of the retarded, advanced and 
causal propagators of all these fields that enable us to define the in and out 
fields of the de Sitter QED (THUS]. 

We must specify that our results based on canonical quantization are dif- 
ferent from those derived by other authors which developed either a de Sitter 
quantum field theory [20] or a quantum mechanics [2T] in axiomatic man- 
ners, trying to avoid the difficulties of the canonical theory by exploiting the 
high symmetry of the de Sitter manifold. Thus, in Refs [21] the tetrad-gauge 
covariance is neglected focusing only on the usual linear representations of 



the SO (1,4) group. However, these representations do not play any role 
in our canonical approach where the the matter fields transform according 
to covariant representations [HI [18]. In other respects, we observe that the 
vector propagators of the canonical theory JT3] are different from the max- 
imally symmetric 5*0(1,4) two-point functions of Refs. [20] even though 
these satisfy the same equations. Moreover, the propagator of the Dirac field 
we derived as a modes sum [T2] was put in a closed form [22] which is no 
longer maximally symmetric. The explanation is that in our approach the 
propagators of the fields with spin get some factors produced by the SL(2, C) 
gauge symmetry which induces the covariant representations of the isometry 
group. 

According to the above arguments we believe that our approach is appro- 
priate for studying QED effects on the de Sitter manifold using perturbations 
as in special relativity [261 [28]. We assume that the QED transitions are mea- 
sured by the same global apparatus which prepares the free quantum states. 
Then the in and out asymptotic fields are free fields minimally coupled to 
gravity but without electromagnetic interaction. Any transition amplitude 
between in and out states has to be calculated by using the scattering op- 
erator S expanded in terms of in fields as in the flat case. However, the 
principal difference is that the de Sitter geometry changes during the infinite 
period of the in — out transitions when the space of this manifold is expand- 
ing from zero to infinity. This could lead to dramatic consequences including 
new quantum effects due to the gravity of the background. 

In this paper we should like to study such simple effects focusing only on 
the non-vanishing in — out amplitudes in the first order of perturbations as 
predicted by the QED on the de Sitter expanding universe. These amplitudes 
are very interesting since they correspond to transitions which in the flat 
limit are forbidden by the mass-shell constraints. On the other hand, we 
know that in QED the non- vanishing first-order contributions are dominant, 
determining the principal properties of these amplitudes. Our main objective 
here is to show that on the de Sitter spacetime there are in — out transitions 
whose first-order amplitudes are non- vanishing and can be expressed in closed 
analytical forms. This enables us to calculate a new dynamical effect of 
particle creation showing that the classical gravitational field can give rise to 
quantum matter during the expansion of the de Sitter universe. Notice that 
the idea of the space expansion generating quantum matter is very old [23] 
but was less studied so far either in the WKB limit [24J or in the quantum 
theory but with a wrong geometric context [25]. For this reason we hope 



that the new results we present here may be of interest. 

We start in the second section presenting a natural version of QED writ- 
ten in moving charts with the conformal time and diagonal gauge. The S 
operator defined as in the flat case is expanded in series of perturbations in 
the momentum representation using the free Dirac and Maxwell fields. In 
the next section the first order in — out transition amplitudes are calculated 
in terms of hypergeometric functions. The section 4 is devoted to the am- 
plitude of the particle creation from the QED vacuum. Concluding remarks 
are presented in the last section. 

2 The de Sitter QED 

The tetrad-gauge covariant quantum field theory on a curved manifold (M, g) 
can be built in any chart {x M } (//, v, ... = 0,1,2,3) where one chooses the 
components of the tetrad fields e^ and e^ defining the local orthogonal frames 
and coframes. These fields are labeled by local indices (with hat) which are 
raised or lowered by the Minkowski metric rj = diag(l, —1,-1,-1) while 
for the natural indices we have to use the metric tensor g^ v = ^ & «e^e^. The 
commutation rules of the vector fields dp = e^d^ define the Cartan coefficients 
as [dp,, dy\ = C^.da- The connection components of any matter field of 
arbitrary spin can be written in local frames using the Cartan coefficients 
and the basis-generators of the finite-dimensional representations of the gauge 
group SL(2, C) transforming the components of the matter fields |S]. 

Assuming that the matter fields of the tetrad-gauge covariant theory are 
minimally coupled to the gravity of M we consider the QED action 

S = I d 4 xy/g (C D + C em + C int ) , g=\deb(g,J)\. (1) 

The Lagrangian density of the free Dirac field ip of mass m has the form [T2] , 



£d = ^ = WD&il> ~ (W) 7 a V] - m^ : , ^ = </> + 7° , (2) 

where the point-independent Dirac matrices 7^ satisfy {7°, 7^} = 2rj a/3 and 
define the basis-generators S a ^ = i[7 a ,7 /3 ]/4 of the spinor representation of 
the SL(2,C) group. The spin connections f A = {{C ao% - C aXo - C oXa )S» X 



give the covariant derivatives of the Dirac field in local frames D& = e^D 



t 1 4A uv\ u,\v v\fi> 



da + fa. The term corresponding to the free electromagnetic field A reads 

i 

C — — — ■ F ~F^ U ■ F - f) A — f) A ("W 

while the fields of the interaction term, 



dm = -e : ip(x)Ye^(x)A u (x)ip(x) : , (4) 

are coupled by the elementary electric charge e corresponding to the U{l) em 
gauge. The notation :: indicates the normal ordering of the products of the 
quantum fields that holds only when the vacuum state is well-defined and 
stable [1]. 

In what follows we consider (M, g) be the de Sitter spacetime where we 
choose the moving chart {£, x} of conformal time, t G (— oo,0), Cartesian 
coordinates and the line element 



ds =- — -r (dt — dx ■ dx) . (5) 

(cut) 2 v / v ' 



The tetrad fields we use have only diagonal components, 



1 g^-rf 1 

"'~' ] uV i i ut 



e° = -cut, e\ = -6\ut, e° = —, e\ = -8) — . (6) 



Notice that here the notation to stands for the Hubble constant of the de 
Sitter spacetime. 

In these frames we have solved the field equations of the free fields mini- 
mally coupled to the de Sitter gravity in momentum representation and we 
separated the positive and negative frequency modes in accordance to the 
choice of the conformal vacuum which is of the Bunch-Davies type [TS] and 
stable on M pp. Therefore, the canonical quantization can be performed as 
in special relativity from which we can take over the well-known methods 
of the QED in flat spacetimes. We note that the de Sitter QED vacuum 
does not represent the state of minimal energy of the whole system since the 
classical gravitational field is still present in this state. Moreover, the classi- 
cal field can change energy with the quantum system without any constraint 
since there the dispersion relation does not work. 

The free Dirac equation can be analytically solved in the gauge ()6]) ob- 
taining the momentum and energy representations with correct orthonor- 
malization and completeness properties [121 [13]. The mode expansion in the 



momentum representation, 

ip{t, x ) = / d?p^2 [ u fA x ) a (p> ° r ) + V^ a {x)a c \p, a) , (7) 

is written in terms of the field operators, a and a c , and the particle and 
antiparticle spinors of this basis, Up >ff and respectively Vp >a , which depend 
on the momentum p and the helicity a = ±1/2. The canonical quantization 
requires the electron (a, a)) and positron (a c , a cj[ ) field operators to satisfy 
the non-vanishing anti-commutators [12] 



{a(p, a), a\p\ a')} = {a c (p, a), a c \p', a')} = S^S^p- p') . (8) 

In the standard representation of the Dirac matrices (with diagonal 7 ) 
the particle and antiparticle spinors read [T2] , 

UfA t, f ) = iN ( .tf ( }:Z7£^^%) ) *• «» 

v*At, *) - <N(. t y- ( -llX^lli ) eH " ■ (10) 

where p = \p\ , H^ 2 ^ are the Hankel functions of indices u± — | ± ip with 
u = — , while 

w =csw?- (11) 

is the normalization constant. The Pauli spinors of the helicity basis, £o-(p) 
and r] a (p) = «cr 2 [£ CT (p)]*, fulfill 

°-p£<t(p) = 2po"& r (p), a-prj a (p) = -2parj (7 (p), (12) 

where a; are the Pauli matrices. The particle spinors have the form 



W^wW-^T). W) 



and satisfy the following properties 



cr -p 



EUp) Up ) + = h,2, E*Z*(?)Z*&) + = — (14) 

(7 CT P 



Similar properties can be deduced for the anti-particle spinors rj a associated 
to £„. 

The free Maxwell field on the de Sitter spacetime can be canonically 
quantized in the Coulomb gauge (with A = and diAi = 0) allowing the 
mode expansion [15] 

A{x) = ( d 3 kJ2 [wz tX (x)a(h,\) + w^ x (x)*a\k,X)] , (15) 

J A 

in momentum representation where the modes functions Wj: A depend on the 

momentum k and helicity A = ±1. The photon field operators (a, cv) satisfy 
the canonical commutation relations 

[a{k, A), a\k', A')] = S X \'S 3 (k - %') . (16) 

The solutions of the free Maxwell equation in momentum representation 
are [U] 

*W*' *) = (2^^ e-^^e x {k) , (17) 

where k = \k\. The polarization vectors E\{k) in the Coulomb gauge must 
be orthogonal to the momentum direction, k-E\(k) = 0, for any polarization 
A = ±1. In general, the polarization vectors have c-number components 
which must satisfy [28J 

e x (k)-e y (k)* = 8 xy , (18) 

Here we consider only the circular polarization and, therefore, we define 
the polarization vectors e±\{k) = -jsi^-e\ — iz-i), in the three-dimensional 

orthogonal local frame {e{\ where k = ke^. 

The Fock space can be constructed now by using the above field operators 
and bearing in mind that the QED vacuum state |0) is a stable one obeying 
a(p,a)\0) = a c (p,a)\0) = a(k, A)|0) = 0. The creation operators give rise to 
the other states of the Fock space according to the standard procedure. For 
example, a state with n\ electrons, n 2 positrons and n 3 photons is denoted 
by 

| ni (p, a, -); n 2 (ff, a\ +); n 3 (k, A)) = J(p, a) ni a c \ff, a') n2 ^(k, A)" 3 |0) . 

(20) 

8 



The vacuum and all the states with different number of particles constitute 
the momentum basis of the Fock space. 

In the presence of the classical gravitational field the reduction formalism 
[26] works as in the flat case if the free fields are quantized canonically |27] . 
This enables us to study transition amplitudes, 

(out; (3....\in; a...) = (in; (3....\S\in; a...) , (21) 

between an in state at t = — oo and an out state at t = 0. We note that 
in the moving chart with the proper time t = — -ln(— ut) these states are 
defined as usual for £ — > — oo and t — > oo respectively. These amplitudes 
have to be calculated using the perturbation method based on the expansion 
of the operator S in terms of in free fields. Recent results indicates that the 
definitions of the in and out fields as well as the mechanisms of the reduction 
formalism and perturbation procedure are similar to those of the flat case 
[27]. For this reason we can use the standard expansion 



S = T exp 



4 J d A x^g CU^inAin) 



(22) 



where the operator products are put in the chronological order with respect 
to the conformal time t. The in fields are free fields which can be identified 
with those defined by Eqs. (J7J and ffT5|) . ip in = ip and A in = A. Replacing 
then in Eq. (BJ the quantities y/g = (out)' 4 and r y^e^(x)A u (x) = —cutj-A(x), 
we arrive at the final formula 

ie Q / : ip(x) 7 • A(x) ip(x) : 



S" = Texp -ie J -^^ : ip(x) 7 • A(x) ip{x) : (23) 

which can be used in applications. 

3 First order transition amplitudes 

The QED transitions which have non-vanishing amplitudes in the first order 
of perturbations involve only three particles, a photon and two Dirac parti- 
cles, which can appear in the in or out states. The allowed transitions are: 
one electron (e~) or positron (e + ) emitting or absorbing one photon (7), one 
photon pair creation or annihilation (7 — > er + e + and e~ + e + — > 7), the 
creation from the QED vacuum (vac) of the triplet e~ + e + +7 and the annihi- 
lation of this triplet into the same vacuum. The transition vac — > e~ +e + + 7 



represents a new dynamical effect of particle creation whose mechanism dif- 
fers from that of the Unruh effect |2J which may be considered as being rather 
of a kinetic origin. The argument is that our effect is due to the QED in- 
teraction and disappears in the limit of the free fields (eo — > 0) where the 
Unruh effect is active. 

The first order amplitudes have simple forms depending only on two 
spinors U and V and the photon mode function w. The photon emission 
e~ — > e~ + 7 has the amplitude 

(p',a',-;k,\\ Si\p, a,-} 

= -ieo J T^vPp'A x ) 7 ' ™k,\( x )* U pA x ) ■ ( 24 ) 

When the photon is emitted by a positron we have to replace Up' iCr , — > Vp j(T 
and Up >ff — > Vp' i<r i. Moreover, if we replace wMwwe obtain the amplitudes 
of the transitions e~+7 —¥ c~ and e + +7 —¥ e + in which a photon is absorbed. 
In the cases of the pair creation, 7 — > e~ + e + , and annihilation, e~ + e + — > 7, 
we find the related amplitudes 



(p,a, -;p ,a',+\Si\ k, A) = -(k, A| S^p, a, -;p , a , +) 

d 4 x — 



ie ° I T^^ U pA x ) 7 • ™k \( x ) v p'y( x ) ■ ( 25 ) 



The dynamical particle creation was studied for the first time in [25] 
where the total amplitude was calculated between an in state at t — > — 00 
and the out state at t — > 00. This means that the conformal time t covers the 
expansion period (— 00, 0] followed by a contraction, for t G [0, 00). In this 
way the contraction cancels the effects due to the expansion, vanishing thus 
the total transition amplitude [25J . In our opinion, the contributions of the 
expansion and contraction periods must be treated separately considering the 
out state at t = in the scenario of the expanding universe. Consequently, 
we find that the transitions vac — > e~ + e + + 7 and e~ + e + + 7 — > vac have 
non- vanishing amplitudes which can be derived by replacing w — > w* in Eq. 

(USD. 

All these amplitudes can be calculated according to Eqs. ©, (TlOl) and 
(TT7|) . After a few manipulation we obtain 

(p',o J ,-;k,X\S 1 \p,a,-) 

10 



.^^JP^S'ip-p'-^CiP^^-exikfUp) 



IQ^/n V k 



x 



r(2,i) 



sign(cr)/V \p,p',-k) + sign(cr') ^- ' {p,p',~k) ,(26) 



(2,1), 



(p,cr, -;p',o" / ,+|5'i| fc, A) 



1dv7t V fc 



X 



r(2,2) 



e*"!^ (p, p', fc) - sign(o-a') e -^P z ' z> (p, p', k) , (27) 



(2,2), 



where we denote 

lt b \p,p\q) 



ds s H^{sp)H^{sp') 



e iqs , a, b = 1, 2 . 



the time integrals of Hankel functions in the new variable s 
obvious properties, 

4 (P, p'a) = I± a) (p',P,q), a, b = 1, 2 , 

i±' 1] (p,p',q) = i^' 2 \p,p',-q), 
i± (p,p,q) = i* {p,p,-q), 



(28) 
-t. The 

(29) 
(30) 
(31) 



indicate that only two types of integrals are independent. Therefore we have 

(2 1) 

nothing to lose if we restrict ourselves to study only the integrals I± ' and 

(2 2) 

I± ' which are involved in the structure of our amplitudes. 

The next step is to evaluate these integrals by expanding them in sums 
of integrals of J-functions as it results from Eqs. (j5Tj) and ( )52|) . Thus we 
obtain 



I± 1] {p,p,q) 



I^' 2 \p,p',q) 



-2 {A±(p,p',q) + C±(p,p',q) 

COSIl 7T/i 



-ie^B ± (p,p / ,q) + ie ± ^B ± (p',p,q)} , (32) 
A {e^A ± (p,p>,q) - C ± (p,p' } q) 

COSh TXjJ, L 



+ie^[B ± {p,p',q) + B ± {p',p,q)}} , (33) 



where the new integrals 
A±(p,p',q) 



dssJ l/± (sp)J u± (sp)e 



/\ Ji s 
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q d 

Q±ip[z - sign(q)iO] , (34) 



7T (j)tf\ 2 dz 

poo 

C±(p,p',q) = / dssJ- u± (sp)J_ u± (sp')c iqs 
Jo 

= ~ ,,3 7 Q^-ifc ~ sign(g)iO] , (35) 

IT (*pp>) 5 cte 

can be calculated straightforwardly using Eq. (l5"3"|) while the integrals with 
indices of opposite signs, 



/■OO 

B±(p,p',q) = / dssJ u± (sp)J- v± (sp')e 
Jo 



,„, _ (lT2i/i)cosh7r/i /p\ 2± ^ 



ttA; 2 (1 + V) U' 



/3 3 1 p 2 p' 2 \ 

x F Mo,l,7±^,^TW /2 , ■ — rwT, ?2 , ■ — rw ,(36) 
\2 2 2 k z + sign{q)tu k z + sign(q)ti) J 



result from Eq. ( 15411 . 

The functions A± and C± are expressed in terms of the Legendre functions 
of the second kind, Q u (z ± iO), depending on the variable 

z = P 2 + P' 2 ~ k \ (37) 

2pp' ' v ; 

which takes values in the domain (—1,1) because of the momentum conser- 
vation in the amplitudes (|26|) and (138]) . Bearing in mind that the Legendre 
functions Q u have a branch cut in this domain we see that the small e which 
assures the convergence of these integrals determines the analytic form the 
Legendre functions given in Appendix B. The functions B± have a more 
complicated structure depending on the Appell hypergeometric functions of 
double arguments F4 [29]. Some technical difficulties could arise here be- 
cause of these functions which are less studied so far. Nevertheless, we have 
all the ingredients we need for calculating the analytical expressions of these 
amplitudes. 

4 The dynamical effect of particle creation 

The examples we analyze now are the related amplitudes of the transitions 

vac — > e~ + e + + 7 and e~ + e + + 7 — > vac that read 

(p, 0, -; V ', o-', +; k, A| Si I 0) = - (0| Si I p, a, -; p ', a', +; k, X}* 

12 



e 



jppi 



5 3 (p + p' + k) Cip) * ■ e x (k)* Va , (p') 



16-y/vr V k 
e^I { +' 2 \p,p', -k) - sign(W) e -^I {2 ' 2 \p,p', -k)] .(38) 



Simple kinetic parameters can be introduced in the orthogonal local frame 
{e*j} where k = —ke^. In this frame we take the electron and positron mo- 
menta in the plane (1, 3) denoting their spherical coordinates as p = (p, a, 0) 
and p = (p',(3,7t) where a, (3 6 (0,7r). Then the momentum conservation 
gives the equations k = p cos a + p' cos (3 and p sin a = p' sin [3 from which we 

deduce 

p sin/3 p' sin a 

k sin(a + (3) 



(39) 
k sin(a + (3) ' fc sin(a + /?) 

Moreover, from Eq. (j3"7|) we obtain z = — cos(a + /3) since the angle between 
p and p is just a + (3. 

The next step is to calculate the matrix elements 



M a , a ,(\)=&(p)a-e x (k)*r ](T ip') 



(40) 



corresponding to our geometry. We use Eqs. (TT31 . the relation r] a = 2cr 2 £* 
and observe that in this case £±i(k)* = 4=(±ei + ie*2) finding the following 



matrices 



M(l) 
M(-l) 



v/2 



- cos f cos 1 
sin | cos | 



- cos -2 sin ^ 
sin | sin | 



sin f sin | — sin | cos f 
cos f sin | — cos f cos | 



(41) 
(42) 



for A = 1 and A = — 1 respectively. 

Finally, we use Eqs. (|34l) . ( J35l) and ( )36l) taking g = — /c for writing down 
the final expression 



(P>, -;p',ct / i +;£,A|S'i|0) 



f 3 o 



-^(p + p' + ^M^A) 



16v/7T fe5 

x {J^(a, /3) + ^(a, /?) - sign(aa') [.F_„(a, /?) + Q-^a, /?)]} , (43) 



where the functions 



J>(a,/3) 



/x(/x 



sin 2 (a + (3) 
2 cosh 2 7T/z sinh 7T/x sin a sin (3 



F [1 — ifi,2 + ifi; 2; cos 



2 a + (3 
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GMiP) 



— sinh TifiF [1 — ifx,2 + if/,; 2; sin 



2 a + (3 



(44) 



2/i — i 



1 



A/sin a sin /3 
7r(l + 4/i 2 ) cosh7r/i sin(a + /3) 

1 



/ sin a \ 2 



+ifi 



\ sin/3 




sin 2 a 



*A*; 



-HO 



sin 2 /3 



sin 2 (a + /3) 'sin 2 (a + /3) 

+ (a i — > /?)] 



-HO 



(45) 



depend only on the angles a and /3 and the parameter /i. 

It is remarkable that the above amplitudes depend on the fermion mass 
and the external gravity only on the parameter /i = — . This parameter 
becomes very small under inflation when the Hubble constant u is extremely 
large. This situation is well approximated by the limit of the amplitude (1471) 
for ji — > 0. Taking into account that in this limit (when v± — > |) the Hankel 
functions are of the form (|50l) we can evaluate the integral 



\im l¥' 2 \p,p', -k) 
which leads to the amplitudes 



li 1 



1 



Ti ypj? k + p + p' — i0 



(46) 



lim(p, a,-;p', a', +; k, A| Si\ 0) 

/J,— >0 



a+/3 



60 S^p + p' + ^^M^iX) C ° S = 



(n-fc) 



cos I cos I 



These amplitudes are non- vanishing only if a' = 
find two non-vanishing amplitudes proportional to 



(47) 



-a. Thus, for A = 1 we 



eo 



2( y 27ikY- 



a + (3 
cos X 



tan ^ 
'-tan i 



for a 
for a 



-a 
-a' 



(48) 



Similar results written for A = — 1 show that all these amplitudes vanishes 
for a = (3 = when e~ and e + have parallel momenta in the same direc- 
tion. However, whether e~ and e + have parallel momenta but in opposite 
directions, i. e. (a = tt, (3 = 0) or (a = 0, (3 = n), we can not use the 
general formula fT4T|) being forced to reconsider the momentum conservation. 
Let us take, for example, the out state with a photon having k = —k e% and 
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A = 1, an electron of parameters p = pe 3 and a = | and a positron with 
p' = (k — p)e 3 provided p > k and er' = — |. Then the resulting amplitude 
for jjl — > is proportional to 

C ° ' (49) 



2(2vr)2pv / A?' 

The conclusion is that under inflation the effect of particle creation is favored 
only when it produces pairs of fermions moving in opposite directions. Hereby 
we can ask if this phenomenon could be one of the mechanisms of separating 
the matter and antimatter between themselves. 



5 Concluding remarks 



We succeeded here to write down the transition amplitudes of the de Sitter 
QED which do not vanish in the first order of perturbations. It is remarkable 
that these transitions are due to the classical gravitational field which changes 
energy with the quantum matter eliminating thus the mass-shell constraints. 
All these transitions represent new quantum effects in external gravitational 
field which were less studied so far [25] . 

The above results were obtained using our gauge- covariant theory of quan- 
tum fields based on canonical quantization which exploits the isometries of 
the curved backgrounds. This allows us to correctly define the spin and to 
derive quantum modes as common systems of eigenfunctions of convenient 
sets of commuting operators. Moreover, in this approach the vacuum state 
is stable such that the canonical quantization, the development of the the- 
ory of interacting fields and the perturbation procedures can be done simply 
as in the flat case. For this reason we believe that our approach represents 
a good tool in calculating quantum transitions on the curved backgrounds 
which have isometries. 

However, in the case of the de Sitter manifold, we did here only one 
step to a long way punctuated by many serious difficulties foreshadowed 
by the analytical forms of our amplitudes which are extremely complicated. 
We can imagine that the next orders of perturbations as well as the renor- 
malization procedures will give rise to new technical difficulties in working 
with special functions and solving complicated integrals. The recent studies 
concerning the regularization of the photon [30J and electron-positron [22] 
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self-energy diagrams in the second order of perturbations confirm this per- 
spective. Therefore, new mathematical methods are needed for solving these 
problems if we want to arrive to a strong theory of interacting fields on the 
de Sitter background, with complete Feynman rules and renormalization in 
any order. 

Appendix A: Integrals of Bessel functions 

The Bessel functions of index | are elementary functions, 

K l (z) = x [^e-\ H?\z) = -uP^e iz , H?\z) = iJ^- e~ iz . (50) 

2 \ ZZ 2 V 7T Z 2 V 7T Z 

The first one helps us to evaluate the integrals (155]) by replacing the expo- 
nential function in Eq. (155]) and using the expansions of the Hankel functions 
in terms of Bessel functions J 



i C2 d 



Hg>(z) = ^ J ^ z )j lJ '^ Z \ (51) 

± COSh(7T/i) 

HJS(z) = f^Mp^M . (52 ) 

± COSh(7T/i) 

We obtain thus two types of integrals which can be put in analytical forms. 
The first integral 

dxxz Ki(cx)J u (ax)J u (bx) =- — 7= - r ~r~Qu-^( u ) j (53) 
2 V27T (06)2 du 2 

depends on the new variable u which obeys 2abu = a 2 + b 2 + c 2 . The second 
integral we consider [35] . 

f°° 3 

/ dxx? Ki (cx)JJax)J-Jbx) 
Jo 2 

simru 5 fa\ v „ /3 „ „ , a 2 6 2 \ . .. 

—^""'(i) i(j.i.i+".i -".-?.-?].(«) 

is solved in terms of Appell hypergeometric functions _F 4 depending on double 
arguments. Both these integrals are convergent for 9ft(c) > 0. Therefore, for 
calculating Eqs. (134]) . (135]) and (136]) we are forced to replace c — > e — iq, 
introducing thus the usual e > which finally tends to zero. 
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Appendix B: Legendre functions 

The Legendre function of the second kind can be written as 

Q v (z±iO) = —?— \e^P u (z) - P„{-z)] (55) 

2 sin 7tz/ L J 

for — 1 < z < 1. The Legendre functions of the first kind, 

P„(z) = f(-u,1 + u;1;?-^) , (56) 



are analytic in this domain being represented by the usual Gauss hypergeo- 
metric function F = 2^1- Hereby we obtain the formula 



d „ , . 711/(1/ + 1) 

az 4 sin txv 



( *'*"F[ I -1/, 2 + 1/; 2;- 

+ /•'(!-!/, 2 + i/; 2; -J-^ 



(57) 



which helps us to calculate the functions A±. (for z/ = ±i/i) and Cj- (taking 
i/ = +z/x — 1). 
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